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Abstract 

Voyager images of the icy satellites of Saturn, Dione and Enceladus, suggest they have been geologically 
active and are not only composed of ice. Recent observations by HST have shown the presence of ozone 
at both Dione and Rhea which also implies the presence of molecular oxygen at these bodies. The Cassini 
Plasma Spectrometer (CAPS) will provide the capability to determine the global composition of these 
bodies by measuring the pickup ions produced by the ionization of their sputter produced atmospheres. 
We will present a model of these atmospheres and associated pickup ions and demonstrate CAPS ability 
to distinguish the freshly produced picked up ions from the ambient plasma. Such ions are expected to 
form a ring distribution that will have a uniquely different energy-angle dependence than the ambient 
plasma ions. In the case of Dione we expect the potential for a moderate strength interaction for which 
both Voyager 1 and Pioneer 11 spacecraft measured ion cyclotron waves centered on the Dione L shell 
and near the equatorial plane. Since Enceladus may be the source of the E-ring, some surprises may be 
encountered during its close encounter with the Cassini spacecraft. In the case of Dione we will show that 
a wake pass at 500 km altitude is more than an order of magnitude better than an upstream pass at 500 
km altitude. Pickup ion detection for minor ion species such as NHy is possible for 500 km altitude wake 
pass but not for 11 500 km altitude upstream pass at closest approach. For navigation reasons a 100 km 
pass is not allowed and therefore it is essential to have a wake pass to maximize the science return for a 
targeted flyby with Dione. The CAPS observations when combined with magnetometer, plasma wave and 
energetic particle observations will allow us to estimate the source of ions into Saturn’s magnetosphere 
due to these two bodies and to characterize the nature of the interaction with Saturn’s magnetosphere. 
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1. Introduction 

In this paper we will be primarily concerned 
with the interaction of Saturn’s magnetosphere 
with its icy satellites Dione and Enceladus. In 
the future we plan to also look at the other icy 
satellites Mimas, Tethys and Rhea; in the case of 
Rhea finite gyro-radii effects will be important. 
The icy satellite interaction with Saturn’s mag- 
netosphere can provide very important informa- 
tion about the bulk properties and surface 
properties of these bodies by detecting the 
pickup ions produced by the ionization of their 
sputtered neutral atmosphere (see Johnson and 
Sittler, 1990). The icy satellites are believed to 
be composed primarily of ice with trace 
amounts of 0 2 , 0 3 and other possible frozen 
volatiles such as CH 4 , NIL CO, , and possible 
dark components such as C and S (see Morrison 
et al., 1984). 

At present, most of our knowledge of the Satur- 
nian system can be traced to the pioneering 
observations of Saturn’s plasma population, 
energetic particle population, magnetic fields, 
plasma waves and remote sensing observations 
from the Pioneer 1 1 and Voyager 1 and 2 en- 
counters. The papers by Frank et al. (1980), 
Lazarus and McNutt (1983), Sittler et al. (1983), 
and Richardson (1986) summarize the first 
observations of Saturn’s plasma environment. 
The energetic particle populations of Saturn’s 
magnetosphere were summarized in papers by 
Fillius and Mcllwain (1980), Vein Allen (1980), 
Fillius et al. (1980), Vogt et al. (1982), Krimigis 
and Armstrong (1982), Schardt and McDonald 
(1983), Krimigis el al. (1983), and Van Allen 
(1984). More recently, Maurice et al. (1996), 
Paranicus et al (1997) and Jurac et al. (2001) 
have reanalyzed the Voyager plasma and ener- 
getic particle data to put it in a more synergistic 
form. Richardson and Sittler (1990) combined 
the Voyager ion and electron plasma observa- 
tions to construct an empirical model of Saturn’s 
plasma environment inside of about 15 R s . In 
conceit Johnson et al. (1989) constructed from 
all these observations a model of a neutral cloud 
torus that surrounds Saturn and is primarily due 


to sputtered neutrals coming from the icy satel- 
lites of Saturn. Richardson et al, (1986) made 
the first attempt at combining all these various 
data sets to construct a semi-empirical model of 
Saturn’s plasma and neutral environment includ- 
ing various collisional and ionization processes 
and radial diffusion which enforced constraints 
on the plasma source strength and the rate of 
radial diffusion. 

This work was followed by the discovery by 
using HST of a relatively large toroidal cloud of 
OH that was surrounding Saturn with densities 
as high as 500 enr 3 near the L shell of Enceladus 
(Shemansky et al., 1993; Hall et al. 1996; and 
Richardson et al. 1998). This discovery put into 
question the original work by Johnson et al. 
(1989) for which it was believed that the icy 
satellites dominated the neutral environment 
around Saturn. In order to reconcile these new 
observations Richardson et al. (1998) invoked 
rapid radial transport with time scales at 6 R s of 
only 5 days while in the original work by 
Richardson et al. (1986) it was 700 days. Results 
by Shi et al. (1995) and Jurac et al. (2001) now 
show that the icy satellites and the observed 
grains in the E ring each contribute less than 
10% of the observed neutral population. The 
detection of ozone, 0 3 , in the surface ice of 
Rhea and Dione (Noll et al. (1997)) for which 
the column densities were ~ 2xl0 16 mol/cm 2 
meant that column densities of 0 2 could be as 
high a 10 19 mol/cm 2 within the icy mantles of 
these bodies. Since the temperature of these 
bodies are ~ 75°K to 100°K (Hand et al., 
1981,1982), there could exist a significant 
gravitationally bound Q„ 0 3 (i.e., since ozone is 
unstable it might not contribute to an atmo- 
sphere around these bodies) atmosphere at 
Dione and Enceladus with the scale heights for 
O, of 1 1 7 km and 227 km, respectively. If the 
O, believed to be trapped in the icy surfaces of 
these bodies, contribute to these atmospheres 
then the peak neutral densities could be as high 
as 8.3xlO n mol/cm 3 for Dione and 4.3xlO n mol/ 
cm 3 for Enceladus. Therefore, a wake pass of 
Dione may see large quantities of pickup 0, + 
and 0 3 + near the limb of Dione; while an up- 
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stream pass at 500 km altitude would probably 
not detect such an atmosphere. In the case of 
Enceladus the atmosphere would be less bound 
with a significant fraction of the 0 2 and 0 3 on 
escape trajectories. We may also expect signifi- 
cant amounts of C0 2 to be present, since CO, 
was detected at the Galilean moons (McCord et 
al,, 1998) and in the atmosphere of Titan 
(Samuelson et al,, 1983) where CO was also 
detected (Lutz et al., 1983). The planned close 
encounters by Cassini of these bodies will allow 
us to develop a more accurate picture of their 
role in neutral cloud production and, thus, in the 
plasma and energetic particle production within 
Saturn’s magnetosphere. 

Since these bodies are deep within the magneto- 
sphere of Saturn and its radiation belts, they are 
expected to have, at minimum, a tenuous atmo- 
sphere due to sputtering induced by energetic 
particle bombardment and bombardment by the 
ambient plasma. The icy satellites should also be 
exposed to micrometeorite bombardment which 
will produce a porous regolith on the their 
surface and the emission of particulates into 
Saturn’s magnetosphere (see Ip, 1997). Knowl- 
edge about the composition of their atmospheres 
can provide information about the surface 
composition (Johnson and Sittler, 1990). In the 
case of Dione the leading hemisphere is brighter 
than the trailing hemisphere (Morrison et al., 
1984) which suggests darkening by plasma and 
energetic particle bombardment. There is also 
evidence of a global tectonic system of fractures 
and faults for Dione (Morrison et al., 1984). In 
the case of Enceladus fractured and smooth 
plains dominate which indicates its surface is 
relatively young ~ 10 8 years (Morrison et al., 
1984). Enceladus is also associated with 
Saturn’s E ring whose brightness peaks at 
Enceladus’ L shell (Baum et al., 1981). The 
evidence for tectonic activity for both bodies 
could indicate the presence of, for example, a 
water-ammonia eutectic which melts at 170°K 
within the interior of these bodies (Morrison et 
al, 1984). The close encounters planned for 
these bodies by Cassini will provide accurate 
measurements of the plasma and energetic 


particle populations at these bodies. This data, 
when combined with the latest laboratory 
measurements of the sputteri ng processes 
(Johnson, 1990; Johnson, 1998; Baragiola et al., 
1999) can provide accurate models of their 
sputtered atmospheres at the time of the encoun- 
ters. 

Saturn’s icy satellites are different than our 
moon since their surfaces are made of ice and 
energetic heavy ions will have sputtering yields 
~ 50 to 1000 H,0 mol/ion (Johnson, 1990). 

z ... 

Because of the large yields, trace refractory 
species are expected to be carried off with the 
water molecules. Although some information 
about the icy satellite composition can be ac- 
quired at large distances, the detected species 
will be partially equilibrated and dissociated 
and, therefore, the origins difficult to determine. 
Furthermore, when sampling the sputtering 
atmosphere of these bodies the freshly picked up 
ions will form unique ring distributions which 
can be readily distinguished from the ambient 
ions which obey bi-Maxwellian distributions. 
Pickup ions formed from what we call ambient 
neutrals will also form ring distributions, but 
they will have accumulated over periods ~10 
hour rotation period of Saturn and are expected 
to have pitch angle scattered to shell distribu- 
tions or bi-Maxwellians before being observed 
(see Vasyliunas and Siscoe, 1976). The atmo- 
spheres of these bodies can be sampled by their 
conversion to ions (photoionization, electron 
impact ionization and charge exchange) in the 
form of pickup ions which can then be measured 
by the ion mass spectrometer (IMS) of the 
Cassini Plasma Spectrometer Experiment 
(CAPS) (see Young et al., 2002 for instalment 
description). The CAPS experiment will also 
provide measurements of the ambient plasma , 
ions and electrons, which will be important for 
estimating the charge exchange and electron 
impact ionization rates for the pickup ions. The 
Magnetospheric Imaging Instrument (MIMI) 
(see Krimigis et al., 2002 for instrument descrip- 
tion) will provide the hot plasma and energetic 
particle observations required to estimate the 
sputtering rates at the time of the targeted 
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flybys. As discussed below the pickup ions are 
expected to form ring distributions which will 
have markedly different energy-angle properties 
than the ambient ions which are typically ap- 
proximated by converted bi-Maxwellian distri- 
butions. We note that the pickup ions with their 
unique energy-angle signature can provide the 
most direct information about the atmosphere 
composition. 

As discussed in Johnson and Sittler (1990) an 
extensively sputtered surface will display the 
following properties: 1) volatile species are 
highly depleted so trace species may be hard to 
detect using reflectance spectra; 2) sputtering 
rates for extensively sputtered surfaces are 
related to bulk composition hence CAPS can be 
more sensitive than reflectance experiment; and 
3) younger surfaces such as that present at 
Enceladus will have higher yields. CAPS, which 
provides a very sensitive measurement of the icy 
satellites bulk atomic and molecular composi- 
tion, is complementary to the information 
provided by the Visible Infrared Mapping 
Spectrometer (VIMS) (see Brown et al, 2000 
for description of instrument) on Cassini which 
will provide spatial maps of surface composi- 
tion. The Composition Infrared Spectrometer 
(CIRS) experiment on Cassini (see Kunde et al., 
2002 for instrument description) will provide 
temperature maps of the icy satellite surfaces 
which will be important for estimating sputter- 
ing rates which are temperature dependent (see 
Johnson, 1990 and Bar-Nun et al, 1985). 

Here we note that Pioneer 1 1 and Voyager 1 
observations show that there may be a strong 
interaction between Dione and Saturn’s mag- 
netosphere which is also supported by the 
simulations presented in this paper. For ex- 
ample, Pioneer 11 (Smith and Tsurutani, 1983 ) 
and Voyager 1 (Barbosa, 1993) magnetometer 
observations showed the presence of ion cyclo- 
tron waves in the vicinity of Dione ’s L shell and 
the frequency of the waves are consistent with 
the dominant pickup ion being a heavy ion such 


as H 2 0 + , OH + or O . These waves, which are 
observed to be confined close to the equatorial 
plane, provide indirect evidence for the forma- 
tion of pickup ions in the vicinity of Dione. 
Galileo results (Russell and Huddleston, 2000; 
Russell et al., 2000) show that peaks in the 
magnetometer spectra give information about 
the composition of the pickup ions being pro- 
duced at a satellite body and will thus comple- 
ment the CAPS observations. If the interaction 
is strong enough then the pickup process will 
also generate an ion beam instability (Ma et al, 
1987) in the plasma. Such instabilities produce a 
spectrum of electric and magnetic waves which 
can be measured by the Radio and Plasma Wave 
Science (RPWS) instrument on Cassini (see 
Gurnett et al., 2000 for instrument description). 
If the satellite has a conductive interior then 
there may be a magnetic signature which could 
be measured by Cassini as was the case for 
Galileo at Europa and and Callisto (see, 

Kivelson et al„ 2000; Zimmer et al., 2000). 

Since the discovery of an internal magnetic field 
at the Jovian satellite Ganymede (Kivelson et 
al., 1 998) we cannot rule out the presence of an 
internal magnetic field which would complicate 
the interaction and infer the presence of a 
dynamo operating in the core of the icy satellite. 
Although we think this unlikely Cassini would 
be able to measure the presence of such a field if 
present and the flyby were close enough to the 
icy satellite. Finally, we note that the satellites 
are not expected to have significant conductive 
ionospheres and that the atmosphere (iono- 
sphere) is optically thin and will not be able to 
prevent the flowing plasma from reaching the 
surface of the icy satellites. The importance of 
these observations during targeted flybys of 
these various bodies by Cassini , emphasizes the 
need for wake encounters as close as possible to 
the icy satellite. This will maximize the science 
return from these targeted flybys which will 
become clearer in later discussions. 
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2, CAPS Instrument Description 
and Capabilities with Respect 
to Targeted Icy Satellite 
Encounters 

The overall capabilities of the CAPS instrumen- 
tation are described in detail by Young et al. 
(2002). The CAPS instrument is composed of an 
Ion Mass Spectrometer (IMS), Electron Spec- 
trometer (ELS) and an Ion Beam Spectrometer 
(IBS). With regard to the targeted icy satellite 
encounters the primary instruments are IMS and 
ELS. Here we focus our discussion on the IMS. 
Papers describing the IMS in detail can be found 
in McComas and Nordholt (1990), McComas et 
al. (1998) and Nordholt et al. (1998). Basically, 
the IMS uses time-of-f light technology with 
cylindrical geometry providing an instantaneous 
2D view of the plasma over an angular range of 
160°; the angular resolution is 8.3°x20°. The 
energy resolution AE/E ~ 17%. The spectrom- 
eter has a medium mass resolution feature (M/ 
AM ~ 8) called straight through, ST, and a high 
mass resolution feature (M/AM ~ 60) using 
linear electric field technique, LEE which can 
also be referred to as an isochronous (time 
focusing) technique. The present estimates of 
the IMS geometric factors are G ~ lO 3 cm 2 -ster- 
eV/eV for ST feature and G ~ K) 4 cm 2 -ster-eV/ 
eV for LEE feature. We will assume these values 
for our simulations. CAPS also has an actuator 
which allows the IMS to have nearly 2 jt cover- 
age in about 3 minutes (l°/sec scan rate). During 
the targeted flybys with an equatorial pass we 
recommend the collimator plane to be in the 
equatorial plane of Saturn and have the actuator 
dither above and below this plane to look for 
pitch angle scattering of the pickup ions which 
are expected to have ring distributions. 

The relatively large geometric factor for CAPS 
allows it to detect ambient and pickup ions with 
high sensitivity. By combining ST and LEF 
measurements and using the mass de- 
convolution capabilities of the IMS Spectrum 
Analyzer Module (SAM) (see Sittler, 1993) the 
IMS will be able to resolve the water group ions 
(ELO + , Off, and ffO + ), organics such as CH 4 + , 


NH 3 + , molecular oxygen 0 2 + , ozone 0 3 + , and 
other species such as MgO + . For example, a 
standard M/Q analyzer would not be able to 
distinguish between 0 + and CH 4 + , but because of 
the non-linear characteristics of the LEF 
technique the CH 4 + ion will breakup in the 
carbon foil and the C + ion exiting the foil will 
appear near the atomic C + ion TOF position but 
shifted in TOF enough to allow it to be resolved 
by the IMS. Our simulations, which are based 
on calibration data, will demonstrate this effect. 

With the collimator in the equatorial plane and 
with an instrument energy-per-charge (E/Q) 
sweep time of 4 seconds (1 eV < E/Q < 42 keV 
with logarithmically spaced 63 E/Q steps), 
CAPS should be able to provide 52 km spatial 
resolution of the pickup ions. Telemetry limita- 
tions confine our measurements of the dominant 
ambient and pickup ions to only 7 ions. But, 
CAPS has the ability to alternate between 7 
different ion species after each energy sweep 
and thus provide the sampling up to 14 ions, 
every 8 seconds. But, in addition to the ion data 
(i.e., 7 to 14 ion species with full energy and 
pitch angle coverage), the IMS will also provide 
TOF versus E/Q spectrograms for both ST and 
LEF data (i.e., data summed over 64 E/Q 
sweeps and collapsed in angle for high sensitiv- 
ity so no angular information is available in this 
data set). This data set will allow us to detect 
both major and minor ions with very high 
sensitivity. The time resolution of this data set is 
about 4 minutes or an equivalent spatial resolu- 
tion of about 3,000 km, which is much larger 
than the dimensions of either Dione and 
Enceladus. Therefore, this data set provides no 
spatial information. We will present simulations 
of both data products in the form of 2D energy- 
angle contours of ion count rate for dominant 
ion species and some minor ion species and 2D 
E/Q versus TOF contours of ion counts for all 
ion species used in our model calculations. The 
ion count rate data is referred to as A cycle data, 
while the E/Q versus TOF data is referred to as 
B cycle data. An A cycle data product is com- 
posed of 8 E/Q sweeps and for the highest data 
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rate there will be 8 A cycles per B cycle. Each A 
cycle is 32 seconds long, while a B cycle is 4 
minutes long, 

3, Model of Sputtered 

Atmospheres for Dione and 
Enceladus 

I n order to construct a model of the sputtered 
atmospheres for either Dione or Enceladus one 
needs the total sputtering source rate and the 
energy spectrum of the sputtered neutrals. For 
simplicity we assume spherical symmetry, but 
before the actual targeted flybys we will con- 
sider models which include asymmetries due to 
ion Mach number, ion and electron pressure 
anisotropies and spatial in-homogeneities of the 
satellite surfaces. The latter can be temperature 
maps, surface properties that the imagers might 
detect and possible compositional spatial varia- 
tions that VIMS might detect . A Monte Carlo 
calculation of the sputtered atmosphere includ- 
ing asymmetries caused by Saturn’s gravita- 
tional field should also be carried out. However, 
to obtain preliminary estimates of the molecular 
species to be detected by CAPS we will use a 
simple spherically symmetric formula given in 
Johnson (1990) (see also Watson, 1982 and 
Sieveka, 1983). For this calculation we could 
use the energy spectrum for the sputtered 
neutrals assuming planar binding (i.e., f(E,cos0) 
= 4lJEcos0/(E+U) 3 ) and spherical binding (i.e., 
f(E,cos0) ----- 2Ucos0/(E+lJ) 2 ), for which we refer 
to Johnson (1990). For our calculations we will 
use the formula assuming spherical binding 
where U is a parameter related to the mean 
energy of the sputtered neutrals. U is also 
referred to as the chemical binding energy of the 
particles in the satellite surface and at present is 
not known with high accuracy. 

For our source term S we will use the maximum 
sputtered flux estimated in Jurac et al. (2001 ) for 
both Dione and Enceladus. In their calculations 
they included the enhancements due to finite 
gyro-radius effects (i.e., effective satellite area) 
and angular enhancement caused by the seccp 
dependence of the sputtered flux where rp is the 


angle of incidence of the impacting ion relative 
to the normal of the satellites surface. Jurac et 
al. (2001) also included the photo sputter flux 
term. Therefore, the sputter satellite source rates 
are S D = 9.6xl0 25 mol/s and S E = 2.1xl(P mol/s 
for Dione and Enceladus, respectively. 

Our sputtered atmospheric model assumes, 
(Jurac et al. (2001)), that the incident ion flux is 
dominated by heavies (i.e., 0 + ) and that the 
dominant sputtered neutral is water molecules. 
We will include additional species in our model 
atmosphere: H 2 , 0 2 , CO, and NH, at the 10% 
level relative to the sputtered flux for water. H, 
and O, are direct sputter products due to ener- 
getic ion bombardment of ice (Johnson, 1990). 
CO, has been seen by the Galileo spacecraft to 
be an atmospheric constituent at Callisto 
(Carlson, 1999) and NH, has been suggested as 
a volatile in the Saturnian satellites because of 
Titan (Stevenson, 1982). Since the ejection of 
trace species is temperature dependent (Johnson, 
1990; Brown et al. (1982)) and depends on the 
satellite composition the choice of 10% is for 
convenience. For example, abundance of the 
minor species such as CH 4 and NH, could be 
lower at the few % level. These results can then 
be scaled based on more detailed knowledge or 
new observations. However, the results given 
here should allow the reader to determine the 
ability of CAPS to detect the minor species for 
either of the targeted flybys. We could also have 
included CH 4 but since the surfaces are oxidiz- 
ing we felt that having CO, was more appropri- 
ate. 

In Figures 1 and 2 we show the result of our 
model calculations for Dione and Enceladus, 
respectively, which shows neutral density 
plotted as a function of height relative to the 
surface. In these calculations we used U = 0.055 
eV for H,0, CO, and NH„ while for H, and O, 
we use U = 0.01 eV. For both bodies the maxi- 
mum water molecule density is about 10 5 mol/ 
cm 3 at their surfaces. In the case of Dione the 
next most important neutral near the surface is 
O,, while far from Dione, in order of impor- 
tance, we have CO,, NH,, and H ,. In the case of 
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Dione Neutral Density 
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Figure 1 . Computed neutral densities of Dione 's sputtered atmosphere as a function of height. Assumed 
species are H,0, H„ 0 2 , CO, and NH ; . The latter four species are minor and assumed to be 10% of water. 
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Figure 2. Computed neutral densities ofEnceladus’ sputtered atmosphere as a function of height. 
Assumed species are 11,0, FI„ O,, CO, and NH,. The latter four species are minor and assumed to be 
10% of water. 



Enceladus, 0, dominates relative to all the other 
minor species for all heights. For the remaining 
minor species in order of importance, we have 
CO,, NH, and H,. To some degree the variation 
of the neutral atmospheric density with height 
can be traced to the assumed energy spectrum of 
the sputtered neutrals and the parameter U. With 
4 second time resolution, CAPS should be able 
to measure the height variation of the neutrals 
from the pickup ion measurements which will 
then allow us to estimate U for the different 
neutral species. Deviations from spherical 
symmetry are expected for radial distances ~ 

Hill Lagrange Sphere which has radius 5.7 R 
(i.e, 3,192 km) for Dione and 3.4 R Enceladus (he, 
850 km for Enceladus). For this paper, we will 
ignore these corrections but when we do the full 
Monte Carlo calculations we will include the 
effects of Saturn’s gravity field. 

4. Estimates of Pickup Ion 
Densities as a function of 
Height for Dione and Enceladus 

4.1 Wake Pass 

In Figure .3 we show an assumed targeted flyby 
geometry where the spacecraft passes through 
the satellite’s co-rotational wake at an unspeci- 
fied height H = (a-r 0 ) at closest approach. The 
reduction of this distance will enhance CAPS 
ability to detect the pickup ions since they will 
have less time to pitch angle scatter and spread 
out in phase space (Wu and Davidson, 1972; 
Hartle and Wu, 1973). Pitch angle scattering 
will tend to transform a ring distribution into a 
shell distribution (Vasyliunas and Siscoe, 1976) 
which will spread the ions in velocity space and 
make the measured fluxes lower at a particular 
energy. Also, pitch angle scattering will tend to 
make the ions move along the magnetic field 
line and thus dilute the pickup ion flux in coor- 
dinate space. The net effect is to reduce our 
ability to measure the pickup ions. An optimal 


height for a wake pass is 500 km since it is 
optimal for the Cassini remote sensing instru- 
ments because of the spacecrafts ability to turn 
in order to remove smearing effects and close 
enough to measure ring distributions by CAPS. 
The figure also shows a fluid element as it 
passes through the satellites atmosphere to the 
observation point where the spacecraft is lo- 
cated. As the fluid element moves through the 
atmosphere it will collect pickup ions as the 
neutral atmosphere is ionized by charge ex- 
change reactions (Table 1), electron impact 
ionization reactions (Table 2) and photo-ioniza- 
tion reactions (Table 3). In order to compute the 
pickup ion densities correctly we will also need 
estimates of the ambient densities and tempera- 
tures of the ions and electrons in the vicinity of 
Dione and Enceladus which are given in Table 
4. Table 4 also gives the ambient densities of the 
various neutral components at Dione and 
Enceladus. All the numbers in Table 4 were 
derived from the papers by Richardson et al. 
(1986) and Richardson et al. (1998). Since our 
model calculations use a fluid approximation we 
are ignoring the finite gyro-radius effects of the 
pickup ions. In the case of Dione the gyro-radius 
for H 2 0 + is about 123 km which is small rela- 
tive to a Dione radius (560 km) and at 
Enceladus the gyro-radius for H,0 + is about 20 
km which is much less than the radius of 
Enceladus (250 km). So, to a good approxima- 
tion the fluid approximation should be accurate. 
Deviations from this will occur on the side 
facing Saturn where the pickup ions can gyrate 
into the bodies surface for heights less than the 
ion gyro-radius. So, for heights near the satel- 
lites surface we would expect to see an asymme- 
try in the observed pickup ion densities. From 
Richardson (1986) the observed plasma velocity 
at Dione is about 46 km/s and at Enceladus is 33 
km/s. The magnetic field strength at Dione is 
about 70 nT while that at Enceladus is 312 nT 
(see Connerney et al, 1983). 
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ICY SATELLITE ELYBY GEOMETRY 



Figure 3. Shows encounter geometry and how we compute pickup ion densities as a fluid element passes through the 
satellites sputtered atmosphere and picks up ions as the neutrals are ionized due to photo-ionization, electron impact 
ionization and charge exchange reactions. 


8 


Table 1. Charge Exchange Rates 


Rates, cm 7s 


Reaction 

Enceladus 

Dione 

Reference 

H + + H H + H + 

1.5E-8 

E5E-8 

Neuman et al. [1982] 

0 + + H 0 + H + 

6.2E-9 

6.3E-9 

Stebbings and Rutherford [1968] 

0 + + OH -» 0 2 + H + 

L2E-10 

E2E-10 

Giguere and Huebner [1978] 

H + + H 2 H + H 2 + 

1.7E-10 

E7E-10 

Tawara [1978] 

h 2 + + h 2 ->h 2 + h 2 + 

6.6E-9 

6.6E-9 

Massay and Gilbody [1974] 

H + + H 2 0 H + H 2 0 + 

5.0E-8 

5.0E-8 

Tawara [1978] 

H + + H 2 0 H + H 2 0 + 

8.2e-9 

8.2e-9 

Ip [1997] 

0 + + H 2 0 0 + H 2 0 + 

2.3E-9 

2.3E-9 

Albritton [1978] 

OH + + H 2 0 OH + H 2 0 + 

L6E-9 

L6E-9 

Huntress [1977], Kim et al. [1974] 

OH + + H 2 -> H + H 2 0 + 

LIE-9 

E1E-9 

Huntress [1977] 

h 2 + + H 2 0 H 2 + H 2 0 + 

3.9E-9 

3.9E-9 

Huntress [1977] 

OH + + OH -> H 2 0 + + 0 

7.0e-10 

7.0e-10 

Ip [1997] 

H 2 + + OH H 2 0 + + H 

7.6e-10 

7.6e-10 

Ip [1997] 

H + + 0 H + 0 + 

5.5E-9 

5.5E-9 

Stebbings et al. [ 1 964] 

0 + + 0 -> 0 + 0 + 

9.9E-9 

9.9E-9 

Stebbings et al. \ 1 964] 

0 2 + + 0 2 0 2 + 0 2 4 

6.4E-9 

6.4E-9 

Banks and Kockarts [1973] 

H 2 + + 0 2 ^ h 2 + 0 2 + 

9.4E-9 

9.4E-9 

Tawara [1978] 

0 + o 2 0 + Oi 

2.1E-9 

2.1E-9 

Albritton et al. [1977] 

OH + + 0 2 -> OH + 0 2 + 

2.0E-10 

2.0E-10 

Bortner et al. [ 1 972] 

H 2 0 + + 0 2 -> H 2 0 + 0 2 + 

5.0E-9 

5.0E-9 

Fehsen feld et al. [ 1 967] 

H + + 0 2 -> H + 0 2 + 

1.17E-9 

E17E-9 

Rudd et al. [1985] 

0 + + C0 2 -> CO + 0 2 + 

LIE-9 

E1E-9 

Giguere and Heubner [1978] 

H + + OH ^ H + OH + 

3.0E-10 

3.0E-10 

Giguere and Heubner [1978] 

H + + OH ^ H + OH + 

2.1e-9 

2.1e-9 

Ip [1997] 

H 2 + + 0 -> H + OH + 

LOE-9 

LOE-9 

Giguere and Heubner [1978] 

0 + + H 2 -> H + OH + 

L6E-9 

1.6E-9 

Huntress \ 1 977], Kim et al. [1974] 

H 2 + + OH OH + + H 2 

7.6e-10 

7.6e-10 

Ip [1997] 

0 + + OH -> 0 + OH + 

3.0E-10 

3.0E-10 

Giguere and Heubner [1978] 

0 + + OH -> 0 + OH + 

3.6e-10 

3.6e-10 

Ip [1997] 

0 + + NH 3 NH 2 + OH 

2.2E-9 

2.2E-9 

Huntress [1977], Oppenheimer [1975]* 

0 + + co 2 -> 0 + co 2 + 

LIE-10 

LIE- 10 

Rudd et al. [1985]* 

H + + C0 2 -> H + 0O 2 + 

LIE-9 

1.1E-9 

Rudd et al. [1985] 

h + + nh 3 -> h+nh 3 + 

5.2E-9 

5.2E-9 

Huntress [1977] 

h 2 + + nh 3 ^h 2 +nh 3 + 

5.7E-9 

5.7E-9 

Huntress [1977] 

OH + + NH 3 OH + NH 3 + 

1.2E-9 

1.2E-9 

Huntress [1977] 

h 2 o + + nh 3 h 2 o + nh 3 + 

2.2E-9 

2.2E-9 

Huntress [1977] 

h 2 + + h 2 o -> h 2 + h 3 o + 

3.4E-9 

3.4E-9 

Huntress [1977] 

OH + + H 2 0 0 + H 3 0 + 

1.3E-9 

1.3E-9 

Huntress [1977] 

H 2 0 + + H 2 H + H 3 0 + 

6.1E-10 

6. IE- 10 

Huntress [1977], Kim et al. [1974] 

H 2 0 + + H 2 0 OH + H 3 0 + 

2.1E-9 

2.1E-9 

Huntress [1977] 


*Estimated 
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Table 2. Electron Impact Ionization Rates 


Rates, cm 7s 


Reaction Enceladus Dione References 


H + e -> H + + 2e 
H + e* -> H + + 2e 
H 2 + e -> H + + H + 2e 
H 2 + e* -> H + + H + 2e 
H 2 0 + e H + + OH + 2e 
H 2 0 + e* -> H + + OH + 2e 
NH 3 + e -» H + + NH 2 + 2e 
NH 3 + e* H + + NH 2 + 2e 
H 2 + e H 2 + + 2e 
H 2 + e* H 2 + + 2e 
NH 3 + e ^ NH + H 2 + + 2e 
NH 3 + e* NH + H 2 + + 2e 
H 2 0 + e H 2 0 + + 2e 
H 2 0 + e* H 2 0 + + 2e 
O + e -» 0 + + 2e 
O + e* O + 2e 
H 2 0 + e -> 0 + + OH + 2e 
H 2 0 + e* 0 + + OH + 2e 
C0 2 + e -> 0 + + CO + 2e 
C0 2 + e* 0 + + CO + 2e 
0 2 + e 0 2 + + 2e 
0 2 + e* 0 2 + + 2e 
OH + e OH + + 2e 
H 2 0 + e -> OH + + 2e 
H 2 0 + e* OH + + 2e 
C0 2 + e -> C0 2 + + 2e 
C0 2 + e* -> C0 2 + + 2e 
C0 2 + e -> CO + + 2e 
C0 2 + e* CO + + 2e 
C0 2 + e C + + 0 2 + 2e 
C0 2 + e* ^ C + + 0 2 + 2e 
NH 3 + e -> NH 3 + + 2e 
NH 3 + e* NH 3 + + 2e 
NH 3 + e NH 2 + + H + 2e 
NH 3 + e* NH 2 + + H + 2e 
NH 3 + e -» NH + + H 2 + 2e 
NH 3 + e* NH + + H 2 + 2e 
NH 3 + e N + + H 2 + H +2e 
NH 3 + e* N + + H 2 +H+2e 

* Hot Electron Component 


1.3E-10 

2.45E-9 

3.1E-8 

3.1E-8 

7.3E-14 

3.3E-11 

1.01E-8 

1.92E-8 

1.66E-12 

7.4E-11 

4.06E-8 

4.35E-8 

2.58E-14 

1.79E-12 

5.28E-10 

5.44E-10 

2.26E-10 

2.2E-9 

5.5E-8 

5.5E-8 

1.35E-13 

1.86E-12 

9.51E-1 1 

9.63E-1 1 

8.6E-11 

8.88E-10 

9.32E-8 

9.3E-8 

2.01E-10 

2.0E-9 

8.2E-8 

8.2E-8 

1.26E-13 

9.72E-12 

1.06E-8 

1.16E-8 

2.17E-12 

7.4E-11 

1.22E-8 

1.45E-8 

2.1E-10 

1.98E-9 

1.5E-7 

1.72E-7 

3.7E-10 

3.16E-9 

1.7E-11 

3.7E-10 

4.61E-8 

4.6E-8 

1.73E-10 

2.36E-9 

1.63E-7 

1.9E-7 

4.45E-12 

7.05E-1 1 

6.04E-9 

6.99E-9 

9.5E-14 

5.92E-12 

3.33E-9 

3.86E-9 

2.66E-10 

2.03E-9 

5.23E-8 

5.26E-8 

3.85E-1 1 

6.6E-10 

4.2E-8 

4.32E-8 

1.71E-13 

9.2 IE- 12 

1.49E-9 

1.52E-9 

1.33E-14 

1.26E-12 

5.44E-10 

5.67E-10 


Lotz [1967] 

Lotz [1967] 

Ip [1997] 

Ip [1997] 

Orient and Srivastava [1987] 
Orient and Srivastava [1987] 
Mark et al. [1977] 

Market al. [1977] 

Kieffer [1969] 

Kieffer [1969] 

Mark et al. [1977] 

Mark et al. [1977] 

Orient and Srivastava [1987] 
Orient and Srivastava [1987] 
Lotz [1967] 

Lotz [1967] 

Orient and Srivastava [1987] 
Orient and Srivastava [1987] 
Orient and Srivastava [1987] 
Orient and Srivastava [1987] 
Banks and Kockarts [1973] 
Banks and Kockarts [1973] 
Richardson et al. [1986] 
Orient and Srivastava [1987] 
Orient and Srivastava [1987] 
Orient and Srivastava [1987] 
Orient and Srivastava [1987] 
Orient and Srivastava [1987] 
Orient and Srivastava [1987] 
Orient and Srivastava [1987] 
Orient and Srivastava [1987] 
Mark et al. [1977] 

Mark et al. [1977] 

Mark et al. [1977] 

Mark et al. [1977] 

Mark et al. [1977] 

Mark et al. [1977] 

Mark et al. [1977] 

Mark et al. [1977] 


30 
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Table 4.Enceladus Torus 


Steady State Densities (cm' J ) 


N ce 

N he 

H H + 

h 2 + o o + 

OH OH + 

h 2 o h 2 o + 

0 2 0 2 + 

100.0 

0.2 

50 25 

0.265* 100 50 

500 10 

100 15 

0.0 0.03* 



Steady State Ion Electron Temperatures (eV) 


T ce 

The 

T(H + ) 

T(H 2 + )* T(0 + ) 

T(OH + )** 

T(H 2 0 + )** 

T(0 2 +)** 

3.0 

100 

11 

22 80 

85 

90 

160 


Dione Tonis 


Steady State Densities (cm"’) 


N ce 

N he 

H H + 

h 2 + o o + 

OH OH + 

h 2 o h 2 o + 

0 2 0 2 + 

25.4 

0.4 

80 4.6 

0.81* 180 11.8 

250 4.45 

30 3.6 

0.0 0.15* 



Steady State Ion Electron Temperatures (eV) 


Tee 

The 

T(H + ) 

T(H 2 + )* T(0 + ) 

T(OH + )** 

T(H 2 0 + )** 

T(0 2 + )** 

5.5 

145 

16.5 

33 105 

111 

118 

210 

* Richardson et al. (1986) 


** Assume temperature proportional to ion mass. 
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The results of our pickup ion densities as a 
function of height for a wake pass is shown in 
Figures 4 a and 4b for Dione and Figure 5a and 
5b for Enceladus. We arbitrarily set the closest 
approach distance equal to the satellite radius 
(i.e., 560 km for Dione and 250 km for 
Enceladus). Referring to Figure 4a for Dione the 
H.,0 + densities reach peak values of about 1 .0 
ions/cm 3 near the surface. FIT has the lowest 
density, while H + and 0 2 + have peak densities of 
about 0.06 ions/cm 3 and 0.08 ions/cm 3 , respec- 
tively. Pickup 0 + is relatively low in abundance 
relative to ambient densities and will probably 
be difficult to detect since the dominant ambient 
ion is expected to be 0 + . The relatively flat 
profile for O' at all heights except near the 
surface indicates we are primarily seeing pickup 
ions produced by ambient neutral oxygen. 

Figure 4b shows pickup ion densities at Dione 
for OFF. COT, NHy, ITO + and NUT. Within 
560 km in height all ions should be detectable, 
OH + dominates relative to the other ions in this 
figure and surprisingly, MX) 4 is relatively abun- 
dant. Near Dione’s surface the total pickup ion 
density is about 1 .4 ions/cm 3 which is about 
5.6% the ambient density and we may expect to 
see a modest interaction at Dione, Here we note 
that because of the detection of ozone, (),. at 


Dione (Noll et al, 1997) we may see significant 
fluxes of pickup Oy and 0, + within a few 
hundred kms above Dione’s surface and there- 
fore a strong interaction similar to that observed 
by Voyager 1 at Titan (Neubauer et al, 1984). 

In the case of Enceladus H,0 + clearly domi- 
nates over all other ions with peak densities near 
the surface of 2.2 ions/cm 3 . The next most 
important ion is Oy with peak densities near 
0.15 ions/cm 3 at the surface. The other ions H + , 
0 + and Hy are of low abundance and will only 
be detectable close to the bodies surface. Again, 
because of the high abundance of ambient H + 
and 0 + , pickup H + and 0 + will be undetectable. 
In Figure 5b the ions IF( ) „ OH + , COy and NHy 
should all be detectable. The ion NHy may be 
detectable close to Enceladus’ surface. Here we 
are using as a rule of thumb for CAPS ability of 
detection of 0.001 ions/cm 3 for ST data and 0.01 
ions/cm 3 for LEF data. This will become more 
apparent when we show our ring distribution 
simulations. 

Figures 4 and 5 also show the height variation of 
the pickup ion densities and as mentioned 
previously this height variation will allow us to 
infer something about the energy distribution of 
the sputtered neutrals. 


Pickup Ion Density for Dione (Wake Pass) 



Figure 4a. Computed ion densities at spacecraft for Dione wake flyby as a function of height for H,0 + , Oy, Hy, 0 + 
and H + ions. The height is equal to the a - r Q parameter in Figure 3. 
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Figure 5b. Computed ion densities at spacecraft for Enceladus wake flyby as a function of height for OH + , 
C0 2 + , NH 3 + , NH 2 + and H.Of The height is equal to the a - r 0 parameter in Figure 3. 


4.2 Upstream Passes 

In Figures 6a and 6b we show the pickup ion 
density as a function of distance Y for an up- 
stream pass with closest approach altitude of 
500 km for Dione. Y is the distance along the 
spacecraft trajectory relative to the center of the 
body (i.e., Y = 0 corresponds to closest 
approach) Jn the case of H 2 0 + we have detect- 
able pickup ion fluxes for the ion measurements 
by CAPS. The H 4 and 0 , + measurements are 
also possible for pickup ions since we can use 
the ST measurements for them. H + and 0 + 
pickup ions will be undetectable because of the 
relatively high ambient densities of these spe- 
cies. Referring to Figure 6b we see that OH + 
(LEF detection), CO + (ST detection) pickup ion 
detection is clearly possible while H 3 0 + (LEF 
detection) detection is very marginal at closest 
approach and probably not possible. If we use 
ST data for OH + and H,(T detection will be 
difficult since ambient 0 + and H ,0 + will domi- 
nate. NH 3 + (LEF detection) is not possible. 


In Figures 7a and 7b we show the pickup ion 
density as a function of distance Y for an up- 
stream pass with closest approach altitude of 
100 km for Dione (i.e,, for navigation reasons 
such a close encounter is not allowed but shown 
to compare upstream pass versus wake pass). 
Here, H,0 + (ST detection), 0 2 + (ST detection), 
and H, + (ST detection) pickup ion detection are 
possible out to several thousand kilometers. H + 
and O pickup ions are not detectable because of 
high ambient densities for these ion species. In 
Figure 7b ions, OH + (LEF detection), and C0 2 + 
(S'F detection) are possible to detect in the ion 
data out to 1 000 km for all. To detect NH 3 + and 
H 3 0 + we would need LEF detection but then 
their peak densities would have to exceed 0.01 
ions/cm 3 . N H e (LEF detection) pickup ion 
detection is not possible. The point should be 
made that we are assuming 10% abundance, 
while the abundances could be a few % or less; 
then detection for many of the minor species 
will probably not be possible. 
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Pickup Ion Density for Dione (Upstream Pass) 



Figure 6a. Computed ion densities at spacecraft for Dione upstream pass as a function of transverse distance 
Y = a (see Figure 3) with 500 km closest approach altitude for H O\ Q, + , FI 2 + , 0 + and H + ions. 


Pickup Ion Density for Dione (Upstream Pass) 



Figure 6b. Computed ion densities at spacecraft for Dione upstream pass as a function of transverse distance 
Y = a (see Figure 3) with 500 km closest approach altitude for OH + , C0 2 + , NH 3 + , NH 2 + and H s O + ions. 
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Pickup Ion Density for Dione (Upstream Pass) 



Figure 7a. Computed ion densities at spacecraft for Dione upstream pass as a function of transverse distance 
Y = a (see Figure 3) with 100 km closest approach altitude for H .O\ OF, H 2 + , 0 + and H + ions. 


Pickup Ion Density for Dione (Upstream Pass) 



Figure 7b. Computed ion densities at spacecraft for Dione upstream pass as a function of transverse 
distance Y= a (see Figure 3) with 100 km closest approach altitude for OIF, C0 2 + , NHF, NH, + and H s O + 
ions. 
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Iii Figures 8a and 8b we show the pickup ion 
density as a function of distance Y for an upstream 
pass with closest approach altitude of 500 km for 
Enceladus. In Figure 8a H,C) + (ST detection), and 
0 2 + (ST detection) pickup ion detection at all 
altitudes are possible. H, + (ST detection) pickup 
ion detection is not possible at all heights. H + and 
0 + are not detectable because of the high ambient 


densities of these species. In Figure 8b, OFF (LEI 7 
detection), H ? 0 + (LEF detection) and NH 3 + (LEF 
detection) of pickup ions are below the detection 
threshold of the instrument. We could try and 
detect OH + and FFO using ST data but will need 
to compete with ambient 0 + and H,0 + . NH 2 + is 
not detectable. 



Figure 8a. Computed ion densities at spacecraft for Enceladus upstream pass as a function of transverse 
distance Y = a (see Figure 3) with 500 km closest approach altitude for ILOy 0 2 + , H, + , 0 + and IF ions. 
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Figure 8b. Computed ion densities at spacecraft for Enceladus upstream pass as a function of transverse 
distance Y = a (see Figure 3) with 500 km closest approach altitude for OFF, C0 2 + , NH 3 \ NFL,* and 
FI 3 0 + ions. 


Pickup Ion Density for Enceladus (Upstream Pass) 
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5 ii Figures 9a and 9b we show the pickup ion den- 
sity as a function of distance Y for an upstream pass 
with closest approach altitude of 100 km for 
Enceladus (i.e., here again a flyby distance of 100 
km is not allowed because of navigation issues but 
is shown to compare upstream passes versus wake 
passes). In Figure 9a H,0 + (ST detection), and C), + 
(ST detection) pickup ion detection is possible at all 
distances. H, + (ST detection) pickup ion detection is 
not possible. The detection of H + and 0 + pickup ion 
detection is not possible because of the large ambi- 


ent ion densities for these species. In Figure 9b. 
OH + (LEF detection) and FFO + (LEE detection) 
pickup ion detection may be possible inside 1000 
km; if we use ST data instead their detection may 
be possible at all heights, but then they have to 
compete with ambient 0 + , OH + and H,0 + . NH, + 
(LEE detection) pickup ion detection is not possible 
at all heights. COT (ST detection) pickup ion 
detection is possible at all heights. NH 2 + (ST or 
LEE detection) pickup ion detection is not possible 
at all heights. 


Pickup Ion Density for Encelodus (Upstream Poss) 



Figure 9a. Computed ion densities at spacecraft for Enceladus upstream pass as a function of transverse distance 
Y = a (see Figure 3) with 100 km closest approach altitude for H,0 + , 0, + , If ', 0 + and IT ions. 


Pickup Ion Density for Enceladus (Upstream Pass) 



Figure 9b. Computed ion densities at spacecraft for Enceladus upstream pass as a function of transverse distance 
Y = a (see Figure 3) with 100 km closest approach altitude for OH + , COT, NH. + , NH, + and H 3 0 + ions. 
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At this point we would like to add the following 
caveat that many of the minor ions for which 
detection using either ST or LET data was ruled 
out as not being detectable referred to the ion 
count measurements or A cycle data, while with 
regard to the B cycle data, which is considerably 
more sensitive, we may be able to detect these 
minor ions. But, such detection would provide 
minimal information about the pitch angle 
distributions of the various ion species which is 
needed to distinguish between ambient ions and 
pickup ions which are expected to be in the form 
of ring distributions. 

5. Ring Distribution Simulations 

Here we derive the expressions used to compute 
the simulated ion count rates for CAPS of the 
pickup ions. These calculations assume that the 
ions are in the form of a ring distribution and 
have not had sufficient time to pitcli angle 
scatter before they are observed by CAPS (see 
Wu and Davison, 1972 and Hartle and Wu, 

1973). In the plasma frame of reference one can 
write the velocity distribution for a fresh pickup 
ion (Hartle and Wu, 1973) 

f (v* ) - — ^ — 5 (v; - v 0L ) 5 (v* - v 0z ) 

2ji v 01 

r the particle velocity in the plasma frame with 
v the mean velocity of the plasma which will be 
in the direction of co-rotation. The component 
perpendicular to the magnetic field is indicated 
by the symbol A and the component along the 
magnetic field is along the z axis. One can then 
re-write the above expression in the appropriate 
form when viewed in the spacecraft frame 
ignoring for now spacecraft motion 


f (v) 


2jt J2. v x V x 


J) (v r - JTv^Vjb (v z ) 


In Figure 10 we show a schematic of the ring 
distribution in the spacecraft frame. We are now 
in a position to incorporate this expression In the 
equation defining the observed count rate for a 
given distribution function f(^). The general 
expression for the count rate that CAPS would 
observe is 

CR = e A jjR(E /Q,6,(j);v) f(v) v • fi v r dv r dv,d<f» 

AvAQ 


for which r(e/q, 6 .f.v) > s the response function of 
the IMS with E/Q the tuned energy-per-charge 
of the IMS, 0 the angle relative to the collimator 
normal of the IMS and (j> the angular sector of 
the IMS. This equation reduces to the following 
expression after substituting the above expres- 
sion for f 


CR = 


e An v, V x A(|> 


In 


Here, e represents the detection efficiency of the 
IMS, A is its effective area of detection, and Atp 
is the angular width of the IMS angular sectors 
which are ~ 20°. We then introduce the geomet- 
ric factor G = 8 AA Q A E / E~ 1 (T 3 cm 2 -ster-eV/ 
eV from which we get 


CR = 


G 

aqae/e 



for which AQ = 0.05 ster, AE/E = 0.17, and 0 < 
v < 2 V . In the limit of v = 2V and setting V 

XX XX ° > 

=50 km/s we get 


m = (vj ]d (v v - pv x V x -v 2 x ) + 6^v y + p P x V x - v x ) 

where V x is the flow speed of the ambient ions 
with x along the co-rotation direction. v is the 
velocity of the pickup ion in the spacecraft 
frame. Finally, this expression can be further 
reformulated into a more useful form 


CR= 1.2x1 if n cps. 

for ST data where n is the density of the pickup 
ion in units of ions/cm 3 . So, it can be seen that 
modest densities can produce very high count 
rates in CAPS if we sample the ring distribution 
function before it can pitch angle scatter and 
spread out in phase space. 
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Figure 10. Ring velocity distribution for pickup ions with drift speed V,. 


Simulations of the pickup ions plus ambient ions 
for various species are shown in Figure 1 1 for 
Dione and Figure 12 for Enceladus (i.e., A cycle 
data) for which we assumed for Dione densities 
of 1.0 ions/cm 3 for H,04 0.08 ions/cm 3 for 0 2 + , 
0.06 ions/cm 3 for C0 2 + , and 0.03 ions/cm 3 for 
NH4 and for Enceladus we assumed densities of 

.p 

2,2 ions/cm 3 for H 4)4 0.15 ions/cm 3 for 0,4 
0.025 ions/cm 3 for CO A. and 0.04 ions/cm 3 for 

z ' 

NH 3 + . These pickup ion densities were derived 
from our wake pass calculations at zero height. 
For the ambient ions we assumed convected 
Maxwellians and used the ion parameters in 
Table 4. For reference we show the energy-angle 
coverage provided by CAPS in Figure 13 a and 


1 3b for Op at Dione and Enceladus, respec- 
tively. It should be clear from Figures 1 1, 12 and 
13 that in most cases CAPS will be able to 
measure the ring distributions of the pickup ions 
for such species as H 2 0 + , 04, CO 4, and NH 3 + 
during a wake pass. This should also be true for 
some of the other ion species such as OH + and 
11.04 As can be seen that the count rates can 
exceed several 1 00 kHz for the pickup ions and 
in many cases should allow us to see them many 
1000 of kilometers above the satellites surface. 
Here, we note that for the molecule NH4 we 
will need to use LEF data, so the displayed 
counts in Figure lid and 1 2d will be an order of 
magnitude lower for LEF data. 
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6. Simulated E/Q Versus TOF 
Spectrograms 

We are now in a position to compute simulated 
E/Q versus TOF contour plots of ion counts for 
both ST and LEF data (i.e., B cycle data) for a 
wake pass only. These simulations will require 
us to integrate along the spacecraft trajectory 
over distances ~ 3000 km. Therefore, this data 
set will not provide information about the radial 
variation of the ion density or information about 
the pitch angle distributions. Considering the 
time and spatial scales of each encounter we will 
get one spectrum before closest approach and 
one spectrum after closest approach. But, the 
advantage of this data set is that it provides the 
energy spectrum of each ion species detected 
and because of its higher sensitivity and ob- 
served TOF spectrum allows us to detect the 
minor ions more easily and detect unexpected 
ion species. The results of these simulations are 
shown in Figure 14a and 14b for Dione and 
Figure 15a and 15b for Enceladus. We display 
the data using the full energy range of the 
spectrometer where we have added nearby 
energy steps as done for the B cycle data (i.e., 

32 E/Q steps). Note, that the spectrograms are 
plotted versus TOF channel number for which 
each channel is about 0.75 nsec wide. We have 
also added Poisson statistical noise to the data 
and accounts for its variable nature. 

Some caveats with regard to the ST simulations 
are required to be said. First, for the H + and H, + 


data we have ignored the importance of the H 
fragment, which can exit the carbon foil; this 
fragment is relatively small when compared to 
the neutral H fragment. At present we do not 
have reliable C + calibration data and are using 
appropriately shifted N + calibration data for our 
C + ions; here we are ignoring the C fragment for 
incident C + ions. For NH + we are also using 
appropriately shifted N + calibration data so that 
the i f fragment is being ignored. For the re- 
mainder of the species 0 + calibration data is 
predominately used and this calibration data was 
contaminated with methane, which was evi- 
dently in the calibration facilities ion source. So, 
this data had C and If fragments. In cases of 
OFF, tif f and H 3 0 + we kept the H fragment 
but deleted the C fragment. For 0 + and 0, + we 
deleted the H‘ and C fragments. For NH, + we 
deleted the C and O' fragments, while for CO, + 
and CO + we deleted the H fragment. At present 
the relative abundance of methane and oxygen 
in the ion source is not known so there is some 
intrinsic error in the relative heights of the 
different fragments for the various species. At 
this time these minor corrections do not impact 
the goals of this paper, but are issues to be dealt 
with when it comes to the final calibration of the 
instrument (i.e., using prototype and engineering 
versions of flight unit) and operation of the 
instrument within Saturn’s magnetosphere and 
surrounding solar wind. 
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Figure 14a. Energy versus time-of-flight (TOF) contour map of simulated ions counts for B cycle ST data 
for Dione wake pass encounter. Vertical lines give approximate peak positions for various incident ions. 
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Figure 14b. Energy versus time-of-flight (TOF) contour map of simulated ions counts for B cycle LEF data 
for Dione encounter. Vertical lines give approximate peak positions for various incident ions. 
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Figure 15a. Energy versus time-of-flight (TOF) contour map of simulated ions counts for B cycle 
ST data for Enceladus encounter. Vertical lines give approximate peak positions for various incident: 
ions. 
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Figure 15b. Energy versus time-of-flight (TOF) contour map of simulated ions counts for B cycle LEF 
data for Enceladus encounter. Vertical lines give approximate peak positions for various incident ions. 
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The simulations of the ST data in Figures 14a 
and 1 5a show that we can resolve the H + and H, + 
in the ST data for which the H + extend to higher 
energy because of its greater mass. The same 
can be said for the water group ions being 
separated from both 0 2 + (i.e., note that in the ST 
data 0 2 + cannot be distinguished uniquely from 
S + ) and C0 2 + , respectively, in the ST data. One 
can also see that the accumulated counts for the 
water group ions can exceed 10 4 counts. For the 
water group ions we see the presence of several 
ghost peaks in the ST data; they occur at chan- 
nels ~ 200, 475 and about 1100. The peaks at 
200-250 channel numbers are probably H ions 
coming out of the carbon foil of the IMS due to 
OH + , H,,0 + , H,0 + and NH, + incident ions which 
have hydrogen in them. The peak at 400 to 500 
channel numbers are probably due to O' ions 
coming out of the carbon foil of the IMS due to 
incident water group ions. The peak at channel 
numbers near 1100 are the ITT ghost peaks due 
to incident water group ions and ammonia ions. 
At Dione H + peaks at around 80 eV and occurs 
at ST channel 140 in TOF with its detectable 
LEF ghost peak at channel 260 (i.e., LEF ghost 
peak caused by H + ion striking LEE stop detec- 
tor and secondary electron from LEF stop 
detector striking ST stop detector). Between 100 
eV to 200 eV the ambient H, + and pickup H 2 + 
combine to form a single peak centered around 
ST channel 195; these ions would also be 
detectable using the LEF data. The water group 
ions at Dione are centered around channel 250 
and 900 eV for the primary peak. 0, + peak 
occurs at around 1 keV and ST channel 850 for 
Dione, while the C0 2 + peak is around 2 keV and 
ST channel 980 for Dione. In the case of 
Enceladus we see a similar result except that 
now the signal is confined below 1 keV, while 
for Dione counts extended up to 3 keV. This can 
be explained by the higher ambient co-rotation 
speeds and temperatures at Dione, and corre- 
spondingly larger pickup energies at Dione. 

In the LEF data for Dione we see a clear peak 
for COT at channel number 700 and E/Q ~ 1-2 
keV. At a slightly lower energy and LEF channel 
number of 950-1000 we see evidence of OT and 


NH 3 + . At around 500 eV we see clear evidence 
of 0 + and H,0 + ions between LEF channels 
1000 and 1100. This data set will allow us to 
measure NH, + , In the case of Enceladus the ion 
peaks are less pronounced, but we should be 
able to measure NHT, FEO + , COT and OT. 

For a more quantitative determination we 
consider several cuts of the ST contour maps at 
specific E/Q. For Dione we consider E/Q = 64 
eV, 724 eV and 1450 eV and are shown in 
Figure 16a. For Encel adus we consider E/Q ----- 64 
eV, 362 eV, and 1024 eV and are displayed in 
Figure 17a. These cuts show the dominant 
presence of H + and HT below 100 eV, while at 
724 eV for Dione and 362 eV for Enceladus the 
spectra are dominated by water group ions and 
finally at 1450 eV for Dione the dominant 
species are 0 2 + and C0 2 + . For Enceladus O T is 
important at 1024 eV and C0 2 + which is impor- 
tant at higher energies is barely seen at 1024 eV. 
At higher energies the peak at 200 to 250 ST 
channels is the H ghost peak due to water group 
ions, at ST channels 400 to 500 we are seeing O 
ghost peaks from the water group ions. Between 
ST channels 1 100 to 1200 we see the LEF ghost 
peaks due to the water group ions. For each 
individual cut of the spectrogram we can apply 
the SAM (Spectrum Analyzer Module) algo- 
rithm (see Sittler, 1993) and de-convolve the 
mass peaks to obtain a more accurate estimate of 
ion counts for each species present in the data. 
The application of this algorithm is no more 
accurate than the accuracy of the calibration 
data. In order to acquire the ion counts for 
various ion species at specific E/Q from this 
data set we need to use the SAM algorithm. For 
this analysis we used a total of 10 ions which 
required us to invert a 10X10 coupling coeffi- 
cient matrix to get a 10X10 inverse coupling 
coefficient matrix (Note, each individual cut 
contains the accumulated counts for 14 ions. For 
our de-convolution algorithm we only used 10 
of those ions and correspondingly only included 
the accumulated counts for those 10 ions.). We 
then bin the data (i.e., sum the ion counts for a 
range of TOF channels centered on ion peak in 
TOE) for 10 ion species and multiply the binned 
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Figure 16a. Simulated ion counts versus TOF for B cycle ST data for Dione encounter. Each curve is a cut 

through the energy-TOF spectrogram at E/Q = 64 eV ( ), 724 eV ( 00 ) and 1450 eV ( ). Vertical 

lines give approximate peak positions for various incident ions. 


Enceladus ST Simulated TOF Spectrum 



TOF Channel 


Figure 17a. Simulated ion counts versus TOF for B cycle ST data for Enceladus encounter. Each curve is a 

cut through the energy-TOF spectrogram at E/Q = 64 eV ( ), 362 eV ( 00 .) and 1024 eV ( ). 

Vertical lines give approximate peak positions for various incident ions. 
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data by the inverse coupling coefficient matrix 
to get the ion counts for each individual ion 
species. The details of applying this algorithm 
are described in Sittler (1993). 

In Tables 5a and 5b we show for Dione and 
Enceladus, respectively, the de-convolved ion 
counts, <C 0UT >, with error bars, a C0UT , for the 
various ion species derived from a subset of the 
ST simulated spectra similar to that shown in 
Figures 16a and 17a. The <C 0UT > is an average 
computed from 40 solutions for which we added 
Poisson statistical noise to the simulated TOF 
spectra. The same can be said for our estimate of 
the standard deviation cr C0UT . The tables also 
show the input counts C IN for ambient ions only 
and ambient ions plus pickup ions. We show the 


ambient only counts to provide a baseline to 
determine our ability to separate the ambient 
ions from the pickup ions. In most cases, 

<C 0UT > « C IN (ambient+ring ) . The results show 
that one can separate the water group ions very 
well and shows the power of our de -convolution 
algorithm which we refer to as the SAM algo- 
rithm. For the water group ions the errors range 
from 1% to 6%. By comparing the ambient only 
and ambient plus pickup ion counts one can see 
that with the B cycle data it will be difficult to 
separate the pickup ions from the ambient ions 
for H + . Hy, and 0 + . For the remaining ions OH + , 
H 2 0 + , H,()\ 0 2 + , COT and C;0 + we maybe able 
to separate them from the ambient ions when 
using the B cycle data when combined with the 


Table 5a. Dione SAM ST Derived ion Counts 


Ion Type 

E/Q (eV) 

Resolution 

C[N 

(ambient) 

ClN 

(ambient +ring) 

< CoUT > 

OCOUT 


1450.0 

ST 

0.0 

0.0 

-1.65 

8.82 

u 

724.0 

44 

0.0 

0.0 

-1.71 

35.6 

u 

64.0 

“ 

2.45E+5 

2.61E+5 

2.61E+5 

637.0 

h 2 + 

1450.0 

ST 

0.0 

0.0 

-1.36 

10.6 

u 

724.0 

44 

66.2 

66.2 

64.8 

36.9 

u 

64.0 

44 

3.91E+4 

4.87E+4 

4.87E+4 

317.0 

0 + 

1450.0 

ST 

2.12E+4 

2.31E+4 

2.33E+4 

3771.0 

44 

724.0 

44 

7.88E+5 

8.11E+5 

8.11E+5 

1.22E+4 

“ 

64.0 

44 

4.48E+4 

4.92E+4 

4.88E+4 

2809.0 

oh + +nhT 

1450.0 

ST 

1.25E+4 

2.44E+4 

2.36E+4 

8618.0 

44 

724.0 

44 

3.06E+5 

3.66E+5 

3.66E+5 

2.66E+4 

44 

64.0 

44 

1.48E+4 

2.86E+4 

2.96E+4 

6096.0 

h 2 o + 

1450.0 

ST 

1.47E+4 

8.90E+4 

8.99E+4 

7728.0 

44 

724.0 

44 

2.51E+5 

6.26E+5 

6.26E+5 

2.29E+4 

44 

64.0 

44 

1.05E+4 

5.86E+4 

5.76E+4 

5154.0 

h 3 o + 

1450.0 

ST 

0.0 

3.83E+3 

3.16E+3 

3.84E+3 

44 

724.0 

44 

0.0 

1.85E+4 

1.83E+4 

1.11E+4 

44 

64.0 

44 

0.0 

2357.0 

2821.0 

2522.0 

0 2 + 

1450.0 

ST 

6373.0 

3.05E+4 

3.05E+4 

286.0 

44 

724.0 

44 

7445.0 

2.17E+4 

2.16E+4 

270.0 

44 

64.0 

44 

120.0 

2250.0 

2247.0 

89.0 

co 2 + 

1450.0 

ST 

0.0 

2.35E+4 

2.35E+4 

293.0 

44 

724.0 

44 

0.0 

7811.0 

7800.0 

184.0 

44 

64.0 

44 

0.0 

1991.0 

1985.0 

71.2 

o 

o 

+ 

1450.0 

ST 

0.0 

681.0 

699.0 

160.0 

44 

724.0 

44 

0.0 

482.0 

482.0 

224.0 

44 

64.0 

44 

0.0 

73.0 

52.4 

79.6 
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Table 5b. Enceladus SAM ST Derived Ion Counts 


Ion Type 

E/Q (eV) 

Resolution 

ClN 

(ambient) 

ClN 

(ambient+ring) 

< CoUT > 

OcOUT 

H + 

1024 

ST 

0.0 

0.0 

-0.22 

5.82 

“ 

362 

“ 

0.0 

0.0 

-6.45 

41.7 

“ 

64 

“ 

2.12E+7 

4.29E+5 

4.29E+5 

869.0 

h 2 + 

1024 

ST 

0.0 

0.0 

3.85 

6.32 

“ 

362 

“ 

456.0 

456.0 

446.0 

48.0 

« 

64 

“ 

3.10E+4 

3.15E+4 

3.15E+4 

279.0 

0 + 

1024 

ST 

4.49E+4 

4.49E+4 

4.51E+4 

3228.0 

“ 

362 

“ 

2.01 E+6 

2.03E+6 

2.03E+6 

2.64E+4 

“ 

64 

“ 

3.64E+5 

3.67E+5 

3.65E+5 

8019.0 

oh + +nh 3 + 

1024 

ST 

1.4E+4 

1.40E+4 

1.36E+4 

7302.0 

“ 

362 

“ 

4.21E+5 

4.51E+5 

4.53E+5 

5.56E+4 

“ 

64 

“ 

6.47E+4 

7.39E+4 

7.94E+4 

1.76E+4 

h 2 o + 

1024 

ST 

2.53E+4 

2.53E+4 

2.53E+4 

6364.0 

66 

362 

(( 

6.24E+5 

1.23 E+6 

1.23E+6 

4.31E+4 

66 

64 

66 

9.56E+4 

2.8E+5 

2.75E+5 

1.53E+4 

h 3 o + 

1024 

ST 

0.0 

0.0 

194.0 

3065.0 

66 

362 

(( 

0.0 

2.03E+4 

2.0E+4 

1.91E+4 

66 

64 

66 

0.0 

6069.0 

8673.0 

7535.0 

o 2 + 

1024 

ST 

0.0 

1.57E+4 

1.57E+4 

177.0 

66 

362 

66 

0.0 

2.3E+4 

2.31E+4 

256.0 

66 

64 

66 

0.0 

5572.0 

5539.0 

122.0 

C0 2 + 

1024 

ST 

0.0 

5565.0 

5545.0 

99.1 

66 

362 

66 

0.0 

3290.0 

3254.0 

183.0 

66 

64 

66 

0.0 

538.0 

536.0 

97.3 

CO + 

1024 

ST 

0.0 

18.5 

-3.8 

83.0 

66 

362 

“ 

0.0 

74.9 

54.2 

295.0 

66 

64 

“ 

0.0 

32.3 

18.8 

152.0 


A cycle data. The primary reason for the ambi- 
guity in the B cycle data is because we add the 8 
angular sectors of the IMS together in the B 
cycle data. In the A cycle data we retain the 
energy-angle dependence of the ion counts. The 
advantage of the B cycle data is its greater 
sensitivity. Well determined ion counts for H + , 
H, + , 0 2 + and CO A are estimated (i.e., errors ~ 


1 %) because these species are well separated 
from each other and the water group ions. NH 3 + 
overlaps with Oil and we are not able to distin- 
guish between the two ions (i.e., one exception 
is the OH + will produce an O' ghost peak and 
NH 3 + will not). The detection of H 3 0 + has a 
large uncertainty ~ 60% due to its co-location 
with the ambient ions 0 + , H,0 + and OH + . 
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The LEF simulated B Cycle data. Figures 1 4b 
and 15b, shows singular peaks for atomic H 4 , 

H 4 H O and ( >'. The simulations also show 
shifted peak positions due to the presence of 
molecules; for example, the figure shows shifted 
peak for C + and 0 + for which the parent mol- 
ecule is COT ; CO + is present but is not re- 
solved because of its low abundance relative to 
the other species. This shift results from the 
molecule breaking up in the carbon foil such 
that the exit speed from the carbon foil of the 
fragments is equal to the entrance speed of the 
parent molecule to the carbon foil. This causes a 
15 ke V CO + incident on carbon foil to have C + 
ion of energy 6,42 keV exiting carbon foil so 
that the fragment has sufficiently low enough 
energy that the detuning effect of the LEF 
section is important (see Nordholt et al., 1998). 
The energy dependence of LEF molecular peaks 
is that the peak position shifts to lower FOP 
with decreasing ion exit energy from the carbon 
foil. In order to calibrate our simulations we 
have taken as a standard, calibration data for N 4 
and combined with the calibration data by 
Nordholt et al. (1998) to extend the coverage 
below 1 5 keV. For E/Q > 15 keV we fit the data 
with a power law in energy (i.e„ TOE = a*E A c + 
b where a,b and c are determined from numeri- 
cal fits to calibration data), while for E/Q < 15 
kV we have fit the data with polynomial of order 
3 (i.e., TOF = ap + bp*E + cp*E A 2 + dp*E A 3 for 
which ap, bp, cp and dp are determined from 
numerical fits to calibration data). For the other 
ions we then multiply the above expression by 
sqrt(M/14) where M is the mass of the fragment 
exiting the carbon foil. In the case of CO + —► C + 
the mass M = 12 and energy E = E co + (12/28) 
which is then substituted into our polynomial 
expression for TOF. In reality more detailed 
calibration data will be required under real 
conditions. These results are accurate enough to 
display the effects expected based on instrument 
response. 

Using the data from Nordholt et al. (1998) 

Figure 6 we have also estimated the probability 
for a C + fragment (P = 0.79) relative to that for 
an 0 + fragment (P = 0.21) for which the parent 


molecule is CO + . A similar calculation can be 
done for COT for which P = 0.65 for C + and 
P = 0.35 for 0 + , Inspection of Figure 14b shows 
no evidence of CO + at TOF ~ 780 which means 
that the peak at about TOF ~ 980 is due to NILT 
One can see from these same figures the shift of 
FIT from H + to lower TOF so that it is well 
separated from the H + peak. We also see a well 
separated peak for NIL, 4 from N + and FLO 4 from 
0 4 . We also see a well resolved peak for (),* — ► 
O 4 ; because of this we can uniquely separate 0 2 + 
from ST While for the ST simulations this 
separation of peaks is not always possible. The 
simulations also underscore the importance of 
having good calibration data. In some instances 
the ST data can better separate the ions of 
different M/Q, while the LEF data is best for 
separating atomic and molecular species of 
similar M/Q. It should also be emphasized that 
in many cases the LEF feature of the IMS will 
provide a unique separation of molecules from 
molecules and atomic ions of the same M/Q. For 
example, the Nordholt et al. (1998) results in 
Figure 6 of that paper show how the CO 4 ~»CT 
and CO 4 -+ O 4 peaks are well separated from the 
N, 4 N 4 peak; for which CO 4 and N, 4 have the 
same M/Q, Also, 40 Ar will be well separated 
from CO , + which has M/Q=44. It should be 
clear from these results that CAPS should be 
able to separate the molecules from the corre- 
sponding atomic species. Again, the main 
difference between Figure 14b and Figure 15b is 
that for Enceladus the ion peaks are confined to 
lower energy relative to that for Dione. Because 
the ions are measured at lower energy for 
Enceladus, the LEF TOF peaks will be shifted to 
lower TOF for Enceladus relative to that for 
Dione. In the case of the ST data the inverse is 
true. Here, we emphasize the fact that the water 
group ions are not well resolved in the ST data 
while for the LEF data they are resolved to a 
greater extent. 

In Figures 16b and 17b we show cuts through 
the LEF E/Q versus TOE spectrograms for 
Dione and Enceladus, respectively, with E/Q = 
64 eV, 724 eV and 1450 eV for Dione and E/Q = 
64 eV, 362 eV and 1024 eV for Enceladus. For 
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Dione we see at E/Q < 100 eV well resolved 
peaks for H + , H, + , NIL" — ► N", H 2 0 + — ► 0 + and 
O". At E/Q = 724 eV for Dione we see well 
resolved peaks for CO," — » C", CO,," —► 0 + , 0 2 + 
-* 0 + , NH," — > N", H,0" — ► 0 + and O. These 
spectra show the power of the LEF feature of 
the IMS in that the species NH ? ", N", H,0" and 
0 + can be separated, while for ST data it would 
not be possible to do very accurately. At E/Q = 
1450 eV for Dione the spectrum is dominated by 
peaks for C0 2 + -+ C + , CO, 1 — » 0 + , 0 2 " —► 0 + , 
NH," — ► N", H 2 0 + — ► 0 + and 0 + . The cuts in E/ 
Q for Enceladus show peaks for H + , H," at 64 
e V and peaks for O," — *■ 0 + , NH 3 " —► N", H 2 0 + 

— ► O' and 0 + . At 362 eV at Enceladus we have 
the same result except the H", It /" ions are 
missing, the COy— » C + peak appears and the 
NH," peak is still well resolved from the 0 2 + 

0 + peak. Finally, at 1024 eV for Enceladus we 
only see peaks due to CO, 1 — ► C", 0 2 + —► 0 + , 

H o+ -+ o+ and O", For Enceladus we do not see 
any obvious evidence of N + in the simulated 
data. 


Similar to what we did for the ST data we are 
now in a position to apply the SAM algorithm 
(Sittler, 1993) to our cuts through the LEF 
spectrograms. In Tables 5b and 5d we show for 
Dione and Enceladus, respectively, the de- 
convolved ion counts with error bars for the 
various ion species derived from the LEF simu- 
lated spectra shown in Figures 16b and 17b. The 
results show that one can separate the water 
group ions fairly well (i.e., uncertainties ~ 1 %) 
with only the H ,( T having relatively large 
uncertainties ~ 15%. Well determined ion counts 
for H + . H , + , 0 2 + and CO, + are estimated because 
these species are well separated from each other 
and the water group ions (i.e., errors ~ 2%-3%). 
For the LEF data we were able to include NH. 4 + 
in our solutions whose peak is well resolved 
from the water group ions; estimated errors are 
~ 3-4%. The data also show that we should be 
able to resolve the pickup ions from the ambient 
ions using the A cycle and B cycle data for ions 
H, + , NHL, OIL, H O", ILCL, 0 2 + , CO;" and 
CO". 



0 200 400 600 800 1000 1200 

TOP Channel 


Figure 16b. Simulated ion counts versus TOF for B cycle LEF data for Dione encounter. Each curve is a 

cut through the energy-TOF spectrogram at E/Q = 64 eV ( ), 724 eV (OO.) and 1450 eV ( ). 

Vertical lines give approximate peak positions for various incident ions. 
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0 200 400 600 800 1000 1200 

TOF Channel 

Figure 17b. Simulated ion counts versus TOF for B cycle LEF data for Enceladus encounter. Each curve 

is a cut through the energy-TOF spectrogram at E/Q = 64 eV ( ), 362 eV (OO..) and 1024 eV 

( ). Vertical lines give approximate peak positions for various incident ions. 
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Table 5c. Dione SAM LEF Derived Ion Counts 


Ion Type 

E/Q (eV) 

Resolution 

ClN 

(ambient) 

ClN 

(ambient+ring) 

< CoUT > 

OCOUT 

H + 


LEF 

0.0 

0.0 

0.1 

0.74 

66 


66 

0.0 

0.0 

0.03 

0.48 

66 


66 

2.45E+4 

2.61E+4 

2.6E+4 

172.0 

h 2 + 


LEF 

0.0 

0.0 

0.23 

1.2 

66 


66 

6.62 

6.62 

6.14 

2.78 

66 


66 

3907.0 

4872.0 

4887.0 

104.0 

nh 3 + 


LEF 

0.0 

271.0 

266.0 

20.7 

66 


66 

0.0 

1360.0 

1355.0 

49.7 

66 


66 

0.0 

313.0 

306.0 

19.3 

0 + 


LEF 

2118.0 

2309.0 

2320.0 

67.3 

66 


66 

7.88E+4 

8.11E+4 

8.11E+4 

313.0 

66 


66 

4475.0 

4923.0 

4914.0 

76.7 

OH + 


LEF 

1246.0 

2168.0 

2148.0 

171.0 

66 


66 

3.06E+4 

3.52E+4 

3.51E+4 

422.0 

66 

64.0 

66 

1476.0 

2542.0 

2578.0 

117.0 

h 2 o + 

1450.0 

LEF 

1470.0 

8899.0 

8878.0 


66 


66 

2.51E+4 

6.26E+4 

6.27E+4 

526.0 

66 


66 

1054.0 

5859.0 

5819.0 

160.0 

H 3 0" 

1450.0 

LEF 

0.0 

368.0 

379.0 

100.0 

66 


66 

0.0 

1851.0 

1873.0 

292.0 

66 


66 

0.0 

236.0 

250.0 

67.3 

o 2 + 

1450.0 

LEF 

637.0 

3049.0 

3039.0 

69.4 

66 

724.0 

66 

745.0 

2168.0 

2163.0 

52.6 

66 

64.0 

66 

12.0 

225.0 

223.0 

22.2 

C0 2 + 


LEF 

0.0 

2350.0 

2337.0 

73.9 

66 


66 

0.0 

781.0 

786.0 

48.1 

66 


66 

0.0 



27.0 

n 

o 

+ 


LEF 

0.0 



30.4 

66 


66 

0.0 

48.2 

46.4 

19.5 

66 


66 

0.0 

7.3 

7.84 

10.0 
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Table 5d. Enceladus SAM LEF Derived Ion Counts 


Ion Type 

E/Q (eV) 

Resolution 

ClN 

(ambient) 

ClN 

(ambient+ring) 

< CoUT > 

CtcOUT 

H + 

1024 

LEF 

0.0 

0.0 

0.0 

0.3 

66 

362 

66 

0.0 

0.0 

0.24 

0.72 

66 

64 

66 

2.12E+6 

4.29E+4 

4.29E+4 

228.0 

h 2 + 

1024 

LEF 

0.0 

0.0 

0.0 

0.44 

66 

362 

66 

45.6 

45.6 

47.2 

9.3 

cc 

64 

66 

3096.0 

3150.0 

3155.0 

70.8 

nh 3 + 

1024 

LEF 

0.0 

0.0 

4.58 

11.4 

cc 

362 

66 

0.0 

879.0 

879.0 

36.8 

cc 

64 

66 

0.0 

265.0 

262.0 

24.7 

0 + 

1024 

LEF 

4486.0 

4486.0 

4470.0 

96.1 

66 

362 

66 

2.01E+5 

2.03E+5 

2.03E+5 

530.0 

66 

64 

66 

3.64E+4 

3.67E+4 

3.67E+4 

233.0 

OH + 


LEF 

1402.0 

1402.0 

1387.0 

114.0 

(( 

362 

66 

4.21E+4 

4.42E+4 

4.41E+4 

627.0 

66 

64 

66 

6471.0 

7125.0 

mmmm 


h 2 o + 

1024 

LEF 

3095.0 

3095.0 

We * 1 

irnw-i 

CC 

362 

66 

6.5E+4 

1.25E+5 

1.25E+5 

763.0 

U 

64 

66 

8679.0 

2.67E+4 

2.67E+4 

325.0 

h 3 o + 

1024 

LEF 

0.0 

0.0 

-6.21 

68.2 

66 

362 

66 

0.0 

2031.0 

1919.0 

470.0 

66 

64 

66 

0.0 

607.0 

594.0 

181.0 

o 2 + 

1024 

LEF 

353.0 

1570.0 

1566.0 

43.1 

66 

362 

66 

538.0 

2304.0 

2311.0 

61.7 

66 

64 

66 

26.7 

557.0 

557.0 

24.9 

co 2 + 

1024 

LEF 

0.0 

556.0 

555.0 

36.5 

U 

362 

66 

0.0 

329.0 

342.0 

26.0 

U 

64 

66 

0.0 

53.8 

55.2 

11.4 

n 

o 

+ 

1024 

LEF 

0.0 

1.85 

5.2 

12.4 

66 

362 

66 

0.0 

7.5 

5.25 

13.6 

66 

64 

66 

0.0 

3.23 

3.96 

6.67 
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7, Summary and Conclusion 

We have presented a model of the sputtered 
atmospheres for Dione and Enceladus where for 
simplicity we assumed spherical symmetry. We 
have used the maximum source rate estimated 
by Jurac et al. (2001) for both these bodies. Here 
the dominant sputtered species is H,0 with a 
10% contribution from (),, H C0 2 and NHL. In 
the case of Dione and Enceladus the maximum 
neutral density at the surface is a few times 10 s 
mol/cm 3 for H,0 with the neutral density vary- 
ing roughly like 1/r 3 with distance from the 
satellite near the body and 1/r 2 far from the 
body. A measurement of the radial variation of 
the neutral density can provide information 
about the energy spectrum of the sputtered 
neutrals (i.e., determine parameter U). By 
measuring the height dependence of the pickup 
ion densities one can reconstruct the height 
dependence of the neutrals. We also speculated 
about the possibility of a gravitationally bound 
atmosphere of 0 3 and 0 2 because of the detec- 
tion of ozone at Dione (Noll et al., 1997). This 
atmosphere, if present, would be considerabOlv 
denser than that estimated for a sputtered atmo- 
sphere and pickup ion densities could be domi- 
nated by 0 2 + and 0 3 + . One could readily detect 
the presence of Argon 40 as radioactive decay 
product which would also be gravitationally 
bound to Dione and Enceladus, respectively. 

We computed the pickup ion densities as a 
function of height by assuming a weak interac- 
tion between the magnetosphere of Saturn and 
the icy satellites atmosphere. We also assumed a 
fluid approximation and neglected finite gyro- 
radius effects which was substantiated by the 
smallness of the ion gyro-radius for pickup 
water ions with respect to the satellite radius. In 
both cases water ions dominated the pickup ion 
density reaching peak values of 1 .0 mol/cm 3 and 
2.2 mol/cm 3 for Dione and Enceladus, respec- 
tively. For a wake pass at the surface of Dione in 
order of most abundant we have OA, OH + , CO A, 
H + , H,0 + , NHA, 0 + , H. + and NHA. In the case of 
Enceladus we have in the order of most abun- 
dant at the surface 0 2 + , H 3 0 + , OH + , NHA and 


COA- At much lower abundance we have NHA, 
HA 0 + and H,A 

We used the pickup ion densities at the surface 
to compute the ring distributions that would 
result and plotted them as a function of E/Q and 
angle in Saturn’s equatorial plane. These simula- 
tions showed the high degree of detection for the 
pickup ions, even when the densities were low. 
This is because the pickup ions obey a delta 
function velocity distribution which translates 
into a high concentration of the ions in phase 
space. The flux peaks at twice the pickup veloc- 
ity. We used these results to show the energy- 
angle coverage provided by CAPS gives us 
ability to resolve the pickup ion distributions 
from the ambient plasma ions. The energy-pitch 
angle dependence of the pickup ions is consider- 
ably different from that of the ambient ions 
which obey converted Maxwellians. 

We also computed simulated energy-TOF 
spectrograms for the Dione flyby and the 
Enceladus flyby. This data is referred to as B 
cycle data (see Young et al. (2000)) and is 
accumulated over several thousand kilometers 
and thus provide no spatial resolution of the 
encounter. The simulations showed that the ST 
data could be used to separate most of the ions 
observed by CAPS but with considerable uncer- 
tainty for some of the minor water group ions, 

N + and NHA- The other ions HA H 2 A OA and 
COA were well resolved by the CAPS ST data. 
While, the LEE data allowed us to resolve the 
water group ions, N + and NH 3 A In the case of 
CO A the LEE data was preferred since it gave a 
unique identification of COA- The same can be 
said for 0 2 + and uniquely separating from SA 
The LEF data also provides the potential for 
uniquely detecting CO + or NA- Therefore, when 
we detect different ion species for high spatial, 
energy and angle resolution we will use a 
combination of ST and LEF data to make the 
ion identification. This will be done by binning 
the ions over specified intervals of TOE and 
then applying the SAM algorithm. For the B 
cycle data we will be able to analyze the energy- 
TOF spectrograms with the maximum resolution 
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in TOF and species and to detect very minor 
pickup species such as FI , + . 

In conclusion, CAPS provides a unique capabil- 
ity in determining the composition of the sur- 
faces of the icy satellites of Saturn (see Johnson 
and Sittler, 1990). Our simulations have shown 
CAPS ability to detect the various ion species 
well above background in most cases, and thus 
infer the bulk composition of the icy satellites. 
Our simulations have not included background 
effects, but is not expected to be important in 
most cases; at most it will just raise the uncer- 
tainty in our determinations. The background 
could be as high as 1 % the peak count rate of 
the dominant species at the energy the dominant 
species is measured. The background wall tend 
to be flat in TOF with an increase toward 
smaller TOF. These corrections, which are also a 
function of the penetrating radiation flux, will be 
included under real measurement conditions. 

For example, if the dominant ion has a peak 
count of 10 5 and the background is 1 % or 10- 
counts, then the uncertainty due to background 
will be ~ 30 counts. (Here we note that this 
background is primarily caused by the scattering 
of ions as they pass through the carbon foil of 
the spectrometer and their accidental detection 
by the start detectors.) With regard to back- 
ground corrections they will be more severe for 
Enceladus than Dione because the radiation belt 
intensities are greater at Enceladus; the back- 


ground is made worse by penetrating radiation 
as demonstrated by the Cassini Earth flyby data. 
Then the background could be more than a few 
percent. At present we are analyzing the back- 
ground corrections due to pentrating radiation 
using the Earth flyby data. We also emphasize, 
because of background corrections, that if we 
want to detect minor ions which could be around 
the few % level relative to 0 + such as NH 3 + will 
be very difficult to detect for an upstream pass 
while possible for a wake pass. In the case of 
NH 3 + we must use LEF data which has a sensi- 
tivity a factor of 10 less than ST data. We have 
also argued that the closer the spacecraft en- 
counters the moon the better, in order to mini- 
mize the diluting effect of scattering which will 
make the ion detection more difficult. A rough 
calculation indicates that pitch angle scattering 
of 5oc will dilute the pickup ions by a factor of 
2; here the pickup ions are confined centrifu- 
gally. We do plan to dither the entrance collima- 
tor above and below Saturn’s equatorial plane to 
look for scattering effects. Finally, our calcula- 
tions with regard to the interaction are not self- 
consistent since we have assumed a weak 
interaction; here we note the possibility for a 
strong strength interaction at Dione, for which 
we note the detection of ion cyclotron waves 
centered on the L shell of Dione by Pioneer 1 1 
and Voyager I . These waves are consistent with 
pickup ions such as the water group ions. 
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